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Abstract 

By means of the nonequilibrium Green function technique, the effect of 
spin-flip scatterings on the spin-dependent electrical transport in ferromagnet- 
insulator-ferromagnet (FM-I-FM) tunnel junctions is investigated. It is shown 
that Julliere's formula for the tunnel conductance must be modified when in- 
cluding the contribution from the spin-flip scatterings. It is found that the 
spin-flip scatterings could lead to an angular shift of the tunnel conductance, 
giving rise to the junction resistance not being the largest when the orienta- 
tions of magnetizations in the two FM electrodes are antiparallel, which may 
offer an alternative explanation for such a phenomenon observed previously 
in experiments in some FM-I-FM junctions. The spin-flip assisted tunneling 
is also observed. 

PACS numbers: 73.40.Gk, 73.40.Rw, 75.70.Cn 
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Spin-dependent electrical transport in magnetic tunnel junctions has received much at- 
tention both theoretically and experimentally in recent years (see e.g. Refs. [Q for review). 
A new field dubbed as spintronics (i.e. spin-based electronics) is emerging. Among oth- 
ers one of the simplest devices in spintronics would be a ferromagnet-insulator-ferromagnet 
(FM-I-FM) structure which is comprised of two ferromagnetic electrodes separated by an 
insulator thin film. In 1975 Julliere made the first observation of spin-polarized electrons 
tunneling through an insulator film into a ferromagnetic metal film, and clearly observed 
14% tunnel magnetoresistance (TMR) for Fe/Ge/Co junctions at 4.2 K 0. In 1995, Mood- 
era et al made a breakthrough that they observed over a 10% TMR for a Co/Al 2 03/Ni 80 Fe 2 o 
junction reproducibly at room temperature ||. Since then, there are a variety of works to- 
ward enhancing TMR in magnetic tunnel junctions, and the TMR > 30% has been obtained 
at room temperature |]J. On the other hand, to understand the spin-polarized tunneling re- 
sults for FM-I-FM junctions, people usually invoke the model based on a classical tunneling 
theory proposed by Julliere 0, in which spins of electrons during tunneling are supposed 
to be conserved, namely, the tunneling of spin-up and spin-down electrons are two quite in- 
dependent processes, and spin-flip scatterings are neglected. Though Julliere's two-current 
model can interpret well some experimental results qualitatively, it still faces to difficulties 
for more complex situations. Actually, in some experiments spin conservation no longer 
holds, and the spin-flip scattering may take effect on the transport properties. There has 
been a number of experiments |l|^] showing that TMR can be very various for different 
barriers, and the inverse TMR can even occur, namely, the resistance when the orientations 
of magnetizations of the two ferromagnets are parallel, is larger than that of antiparallel 
orientations. It appears that the spin-flip scatterings might not be ignored in these situa- 
tions. Recently, Vedyayev et al investigated a model including impurities in the middle of 
the barrier, and considered both cases of spin conserving scattering and spin-flip scattering 
0. Besides, Jaffres et al measured the angular dependence of the TMR for transition-metal 
based junctions, and observed that the angular response is beyond the simple cosine shape 

In this paper, based on a microscopic model and using the nonequilibrium Green function 
technique we shall give a more general expression of the angular dependence of the TMR 
for FM-I-FM junctions by including the effect of spin-flip scatterings. It is found that the 
effect of spin-flip scatterings gives rise to a correction to the formula of the usual tunnel 
conductance. Specifically, we have found that the spin-flip scattering induces a phase shift, 
leading to the tunnel conductance not to be the smallest when the magnetizations of the 
two FM electrodes are antiparallel, which may provide an alternative explanation for the 
previously experimental observation of the angular shift in some FM-I-FM tunnel junctions. 
In addition, it has been shown that the spin-flip scattering could also lead to the inverse 
TMR effect under certain conditions though it is not the only factor, and the spin-flip 
assisted tunneling is also observed. 

Let us consider a magnetic tunnel junction consisting of two ferromagnetic films separated 
by an insulator thin film. A steady bias voltage is applied to the junction. The relative 
orientation of magnetizations in the two ferromagnets is characterized by the angle 9. The 
Hamiltonian of the system reads 

H = H L + H R + H T} (1) 
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with 

ka 

H R = - aM 2 cos 9)c\ a c qa - M 2 sin flc^c^] , 

qa 

H T = l T kq a ko C q°' + T kq * C \a' a ka], 

kqaa' 

where a ka and c ka are annihilation operators of electrons with momentum k and spin a 
(= ±1) in the left and right ferromagnets, respectively, = e L (\s)—eV—aM 1: Mi = smbL.^ 
M 2 = 9flB 2 hR , g is Lande factor, yU B is Bohr magneton, hi(R) is the molecular field of the 
left (right) ferromagnet, e L ( R )(k) is the single-particle dispersion of the left (right) FM 
electrode, V is the applied bias voltage, T%? denotes the spin and momentum dependent 
tunneling amplitude through the insulating barrier. Note that the spin-flip scattering is 
included in when a' — a — —a. It is this term that violates the spin conservation in 
the tunneling process. By performing the u — v transformation, c qa = cos |6 gCT — a sin §6^, 
c\ a = cos \b\ a - asm §6^, H R becomes H R = J2 q a £qab\ a b qa , with = e R (q) - aM 2 , and 
H T becomes H T = Y^kqaa 1 T kq {cos § a\j) qa i — a' sin ^a\j) q ^') + h.c. The tunneling current 
has the form of 

hiY) = e(N L ^ = -|Re^Tr CT {lVG< (*, *)}, (2) 

/ T TT rpU \ ( o _ sin e \ 

where Q kq = T kq • R with T kq = ft ft and R = } e * , and Tr a stands 

\ L kq L kq ) \ bU1 2 LOb 2 / 

for the trace of matrix taking over the spin space. The lesser Green function, G kq (t, t'), is 
defined in a steady state as 

-f) -f)" 



^-(^(t-to Gi!'<(t-t')J' (3) 

with G™ q ^(t — t') = i (a\ a (t')b q(7 r(t)^ . To obtain the lesser Green function, one needs to 
introduce a time-ordered Green function G* fc as 

6 w-o=($!::g;$j!:2). « 

with — i') = — i ^T{& go -(t)a^ CT , (t')}) • By using the equation of motion, we get 

G;*(e) = E F U e ) n Uk( £ ). (5) 
?' 

with 

F^(e) = £ R (e) 5 qq/ + £ G^O^g^e), (6) 
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where the use has been made of the Fourier transform of the time-ordered Green function 
Ggfc(^) = / dte 16 ^'^ G gk (t — t'), and g,l L (e), g\i R {£) are the time-ordered Green function of 
the left and right FM electrodes for the uncoupled system, respectively. By applying the 
Langreth theorem || to Eq.(||), we get 

G< q (e) = £[F^,( £ )nt 9 ,g<,( £ ) + F<,(e)n{ ql g a kL (e)}, (7) 

q' 

where the superscript "r(a)" denotes the "retarded (advanced )" Green function. The tun- 
neling current becomes 



h{V) = ^Re(E{^(^[F^( £ )fit g ,g< L ( £ ) +F< qt (e)nt ql gl L (e)])}). (8) 

k,q,q' 



To get useful analytical result we may assume for simplicity that the tunneling amplitude 
Tkq is independent of the momentum like the conventional consideration |]|, but depends 
on spin. This suggests that apart from inclusion of the spin-flip scattering we have supposed 
that the tunneling amplitude of electrons for the spin-up channel differs from that of the 

spin-down channel. As a result, T kq becomes T = ( _, _, ). In addition, the elements of T 

are assumed to be real. Up to the first-order approximation to the Green function F qq/ (e), 
we get 

h{V) = |Re j ^[/(e)-/(e + eV)]Tr ff Fr e// ( e ,V)], (9) 

where /(e) is the Fermi function, 

T e// (e, V) = 2tt 2 T • R • D R (e) • R f ■ • D L ( £ + eV), (10) 

and Dl(_r) (e) is a 2 x 2 diagonal matrix with two nonzero elements being the corresponding 
density of states (DOS) of electrons with spin up and down in the left (right) ferromagnet. 
For a small bias voltage V, we obtain the tunneling conductance 

2e 2 

G = X T e// , (11) 

where T e ff = Tr a \Re{T e ff{ep, V = 0))], and Bp is the Fermi energy. In comparison to the 
Landauer-Biittiker formula, one may find that T e ff can be regarded as an effective tunneling 
transmission coefficient which includes the contribution from spin-flip scatterings. Eq. (|TT|) 
can be explicitly rewritten as 



G = G [l + P 2 JP? + Pi cos(d - B f )l (12) 



where G = ^[(Tf + Tf)£> iT + (Tf + Tl)D Li ]{D m + D m ), P x = i^gjg^g^g^ , 

p - = SsSt> ^ = ^wSm^fk-^ tand f = ft' D ^ = D ^ e + Ml + eV ^ ^ = 

D L (e-M 1 + eV), D m = D R (e + M 2 ), D Rl = D R (e-M 2 ), and D L{R) is the DOS of electrons 
in the left (right) ferromagnet. One may observe that there is an angular shift induced by 
the spin-flip scatterings, as to be discussed below. 
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Now let us look at the angular dependence of the conductance G. When the spin-flip 
scattering is neglected, i.e. T 2 = T 3 = 0, and if we further assume 7\ = T 4 , we recover 
the conventional expression for the conductance G = Gq(1 + P1P2COS8), which is familiar 
in literature 0,0], where G = ^T 1 2 ( J D L? + D Ll )(D R1 + D Rl ), and P 1 = g^;gg is the 
usual polarization of the left ferromagnet, and P3 = 0. When T2 = T3 = but T± 7^ T4, 
which implies that even if the spin-flip scattering is ignored, but the tunneling amplitude of 
electrons for the spin-up channel is different from that for the spin-down channel, we can 
also get in this situation an expression 

G = G?o(l + P[P 2 cos 9), (13) 

where G' = ?£{T?D L1 + TlD L[ ){D m + D m ) and P[ = ^g^lg^j . Although it looks 
seemingly like the conventional form, it is clear that the difference of the tunneling ampli- 
tudes for the two independent spin channels can still alter the magnitude of the conductance 
and the polarization as well. 

On the other hand, the angular dependence of the tunnel conductance without consid- 
ering the spin-flip effects, as mentioned before, is well known: 

9 9 

G = G P cos 2 - + Gap sin 2 -, (14) 

with Gp the conductance for parallel orientation of magnetizations in the two FM electrodes, 
and Gap the conductance for the antiparallel orientation (see e.g. Refs. [p]J|). While in the 
present case, namely, with inclusion of the effect of spin-flip scatterings we find from Eq.(12) 
that the angular dependence of the conductance becomes 

G = Gicos 2 ^ + G 2 sin 2 | + G 3 sin0, (15) 

where G x = ?£{D m {T?D L1 + TiD Ll ] + Dp^D^ + T 2 D Ll }}, G 2 = ^{D* T [Tf L> iT + 
TlD Ll ) + D Rl [TfD Lr + TiD Ll }}, and G 3 = ^{D R ^ - Dr^T^D^ + T 3 T 4 D Ll ). One may 
see that apart from the conventional cos 2 ~ and sin 2 | terms there is an additional third 
term proportional to sin 9. Here we should point out that G\ is the conductance in the 
case of the parallel alignment (9 = 0) for the magnetizations of the two ferromagnets, G 2 
is corresponding to the antiparallel case (9 = ir), and G3 gives an additional term for the 
noncollinear case (9 7^ and rr), which disappears in the collinear cases. Certainly, these 
three coefficients Gi (i = 1, 2, 3) contain the contributions from the spin-flip scatterings 
characterized by T 2 and T 3 . It is the effect of spin-flip scatterings that enables the tunnel 
conductance not to be at the minimum when the magnetizations of the two ferromagnets 
are antiparallel. This is understandable, because the spin-flip scattering process violates 
the spin conservation and can enable electrons in the spin-up band of one FM electrode 
tunneling through the insulator barrier into the spin-down band of another FM electrode, 
and vice versa, thereby giving rise to a phase shift, as shown in Eq. (|T2|) . It is emphasized 
that this shift will disappear when the effect of spin-flip scattering is neglected. Therefore, 
Eqs. (|T2|) and ( [To] ) can be viewd as a generalization of the conventional expression for the 
tunnel conductance [see Eq.(|T4|)l. It is interesting to note that the phenomenon of such 
an angular shift has been experimentally observed for a CoFe/Al 2 03/Co tunnel junction, 
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as presented in Ref. || (see Fig. 4 therein), where the maximum of the junction resistance 
appears at = 200°, not 180°, implying the angular shift Of — 20°. Although the authors 
of Ref. m did not mention the reasons why such an angular shift occurs in this FM-I-FM 
junction, in accordance with the aforementioned analysis we may attribute this phenomenon 
possibly to the effect of the spin-flip scatterings. If this is acceptable, we can in turn infer the 
magnitude of the effect of spin-flip scatterings. To show it explicitly, let us assume T\ w T4 
and T2 ~ T3 for simplicity, and define a parameter 

7 = |, (16) 

which characterizes the magnitude of the effect of spin-flip scatterings. The angular shift Of 
versus the parameter 7 is plotted in Fig. 1. It is seen that Of is monotonously increasing with 
increasing 7. When 7 approaches to 1, Of = 90° which can be obtained from the expression 
of Pi (see below). If 7 > 1, then P 1 can be negative, leading to Of larger than 90°. To 
see more clearly the effect of spin-flip scatterings on the conductance, we note that Go i n 
Eq. (|i~2|) can be written as G = (1 + 7 2 )G . In this case, Pi = j^sPi, and P 3 = j+s- 
The 7-dependence of the conductance is presented in Fig. 2 (a). It can be seen that the 
conductance decreases with increasing 7 when the magnetizations in the two ferromagnets 
are parallel, while it increases for the antiparallel alignment. As 7 > 1, one may observe 
that G{0 = 7r) > G{0 = 0), suggesting that the inverse TMR may occur. For the case of 
noncollinear alignments, with increasing 7 the conductance first increases rapidly and then 
decreases, and some peaks appear around 7 « 1. 

We come to consider the effect of spin-flip scatterings on the TMR. Recently, a number 
of experiments [|j for a few FM-I-FM tunnel junctions show that if the insulator thin film 
is the material which differs from A1 2 3 , such as SrTi0 3 , Ceo.69Lao.31O1.s45 and so on, the 
TMR, defined as usual as 1 — G(0 = ir)/G(0 = 0), will be negative under certain conditions, 
which means G(0 = 71) > G(0 = 0), exhibiting the so-called inverse TMR effect. It is 
generally believed that this effect may originate from the electronic states at the interface 
between a ferromagnetic layer and an insulating layer |4j] which could give rise to the density 
of states in the minority spin band larger than that in the majority spin band at the Fermi 
level. However, one may see below that the spin-flip scattering can also contribute to the 
inverse TMR. From ( |T2"D we find that the TMR still has the apparently standard form 

9 P P 

TMR = 12 , (17) 
I + P1P2 1 1 

but Pi, defined after Eq. (0) and containing the contribution from spin-flip scatterings, 
differs from the conventional polarization Pi. In Ref. 0, Moodera et al calculated the TMR 
for the CoFe/A^Os/Co junction according to the Julliere's formula. The calculated result is 
27%, while the experiment value is 24% at 4.2 K. If we adopt this experimental data, we can 
calculate the contribution of the spin-flip scatterings which is characterized by the parameter 
7. The obtained result for 7 is about 0.28, where the spin polarization of electrons is taken 
as 47% for CoFe and 34% for Co ||. It shows that the spin-flip scatterings might have a 
considerable effect on the electrical transport of this tunnel junction. On the other hand, 
if we take the contribution from the spin-flip scatterings into account, then we can apply 
our formula to estimate the value of TMR, which could be closer to the experimental result 
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than using Julliere's formula. If the tunneling amplitudes satisfy a certain condition, P\ can 
be negative, depending on the difference between (T-j 2 — T%)D L ^ and (T 4 2 — T%)D L ^ resulting 
in that the TMR can be negative. The 7-dependence of the TMR is depicted in Fig. 2 (b) . 
It can be seen that the TMR decreases with increasing 7, and becomes negative for 7 > 1. 
This can be understood in the following. For < 7 < 1, the tunneling amplitude for the two 
independent channels (T\) is larger than the tunneling amplitude for the spin- flip channel 
(T2). Owing to the spin-flip scatterings the polarization Pi becomes effective and small, 
leading to decreasing of the TMR. When 7 > 1, i.e., the tunneling amplitude for the two 
independent channels is smaller than that for the spin-flip channel, the spin-flip scattering 
dominates in the tunneling process, implying that the electrons with spin up in the left FM 
electrode can tunnel through the insulator barrier to occupy the states of electrons with 
spin down in the right FM electrode via the spin-flip mechanism, thereby giving rise to 
contribution to the inverse TMR effect. 

The angular dependence of the TMR can be understood from the following definition 



TMm = -> = ™ + W - "A . (18) 

When 9 = Of, the TMR goes to its maximum which will be denoted by TMR{6f) here- 
after. In Fig. 3, the 7-dependence of the maximum TMR, i.e. TMR(6f), is presented. A 
remarkable property is that TMR(6j) has a peak at 7 ~ 0.82, which migh be a result of 
the spin-flip assisted tunneling. It is seen that TMR(9f) approaches to a constant when 7 
increases to a large value. No matter how large 7 is, TMR(9f) is always positive. 

The above discussion is based on the result up to the first-order approximation for the 
Green functions. To include contributions from higher-order Green functions, it is better to 
consider them in a Keldysh space. We introduce the nonequilibrium Green function in the 
Keldysh space as || 

^^-(SSoSfto)- (19) 

where G* fc (£, f) is the time-ordered Green function defined as before, G^ ,(t, t') is lesser 
(greater) Green function, and G* fc (t, t') is the antitime-ordered Green function. The elements 
of G * fc (t, f) and G> k (t, f) are given by 

Gg% = -i<f m*)«UO]>, (20) 
Gtf'^t, = -i(M*)aUO)- 

According to Eqs. (^) and ([]), we can make a sum of Feynman diagrams shown in Fig. 4. 
The result is 

G kq (e)=g R (e)£(e)g L (e), (21) 

where — iS(e) = — i&r 3 (l — rj)' 1 , fj = gL(£)^%gR(^)^%, 17 = T • R, and f 3 is a Pauli 
matrix. From Eq. (fj), the current can be rewritten as 
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LP f (IF 

h(V) = - - Re J -Tr a {(nG qk (e)) 12 }. (22) 

After a tedious calculation, we get 

p r (If 

I L {V) = -j-Rej — rr{F 2 (f 3 -if 2 )F 1 T e// (e i y)[(l+T e// (e,V)/2)7b 

+(/(e + eV) - /(e))(f 3 + if 2 )T e// (e, V)]- 1 }, (23) 

where Tr denotes the trace over the spin space and the Keldysh space, f is the unit matrix, 
fj (i = 1,2,3) are Pauli matrices, and 



« _ , 1-2/(5) V « _ /^/(e + eV) 
' 1 2/(e)J'* a I i_2/(e + eV) 



Eq. (|23]) can be rewritten in a compact form 
4e 



4e /" 

I L {V) = —ReJ de[f(e + eV) - /(e)]Tr ff {[A(e, V) - 1] ■ A(e, V)}, (24) 
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where A(e, 7) = • (1 + ^1)-^ and T e// (s, ^) = ( ^ ^ ^ 

matrices in spin space, and = 1,2,3,4) are the elements of the effective transmission 

matrix T e ff(e, V) given in Eq. (|i~0|) . In principle, Eq. (P4|) gives the current including more 
corrections from the spin-flip scatterings, and the tunnel conductance can be obtained by 
G = dIi(V) / dV . For a small bias voltage, we get 

G = ^f eff , (25) 

, [8T° +(Ti ( £F )+T 4 ( £F ))(2+T0 )] 

where Te // p+^Z+V^^)^^))]" with 

T e °jj = 27r 4 D i |(ep)Z) i |(e i ?)Z)j ? |(e^)£) jR |(e:^)(T 1 T4 — T 2 T 3 ) 2 , can be viewed as the effective 
transmission coefficient. Compared to Eq. (|Tl"D, Eq. ( p5|) includes more corrections from 
the spin-flip scatterings. Although the angular dependence of G{6) determined by Eq. (|25|) 
looks more complex in form than one presented in Eq. (|12|) , the behavior of G(8) versus 9 
is found to be qualitatively similar to those shown in Fig. 2 (a) for a given 7. 

In summary, we have investigated the effect of spin-flip scatterings on the spin-dependent 
electrical transport in ferromagnet-insulator-ferromagnet tunnel junctions by using the 
nonequilibrium Green function technique. When the effect of the spin-flip scatterings is 
taken into account, the frequently used Julliere's formula for the tunnel conductance must 
be modified, though the form looks seemingly similar. It is found that the spin-flip scat- 
terings could lead to an angular shift of the tunnel conductance, giving rise to the junction 
resistance not being the largest when the orientations of magnetizations in the two FM 
electrodes are antiparallel, which is quite consistent with the experimental observation in 
CoFe/A^Os/Co tunnel junctions. As the Julliere's formula overestimates the value of the 
TMR, our derived formula with inclusion of the effect of spin-flip scatterings could esti- 
mate the TMR value closer to the experimental result, as discussed above. It is found that 
the spin-flip scattering could also lead to the inverse TMR effect under certain conditions, 
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though it is not the only factor. The phenomenon of spin-flip assisted tunneling is clearly 
observed. When including high-order terms of the Green function, the angular dependence 
of the tunnel conductance is qualitatively similar, although the form looks more complex. 
Finally, we would like to mention that our present derivation can be readily extended to 
other magnetic junctions. 
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FIGURE CAPTIONS 

Fig. 1 The angular shift 9f versus the parameter 7, where the mass of single electron 
in both ferromagnets are assumed as unity, the molecular fields in the two ferromagnets are 
supposed to be the same and taken as 0.7eV, e F = 1.5eV, and the coupling parameter T x is 
chosen as 0.01 eV. 

Fig. 2 The 7-dependence of the tunnel conductance (a) and the TMR (b), where the 
parameters are taken the same as those in Fig. 1. 

Fig. 3 The 7-dependence of the maximum tunnel magnetoresistance TMR(6f), where 
the parameters are taken the same as those in Fig. I. 

Fig. 4 Feynman diagrams for the Green function G qk (e). 
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Fig.2 Zhu et al 
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